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The main subject of these lectures is the Nonperturbative Wave Functions. We 
describe some nonperturbative methods (like QCD sum rules, dispersion relations, 
duality etc) in order to study this object. We also consider some application of 
the obtained results, such as form factors, inclusive amplitudes and diffractive 
electroproduction. Finally, we discuss the instanton liquid model which may help 
us to understand the success of the constituent quark model. 



1 Introduction 

The problem of bound states in the relativistic quantum field theory with large 
coupling constant is, in general, an extremely difficult problem. Understand- 
ing the structure of the bound state is a very ambitious goal which assumes 
the solution of a whole spectrum of tightly connected problems, such as con- 
finement, chiral symmetry breaking phenomenon, and many others which are 
greatly important in the low energy region. Fortunately, a great deal of infor- 
mation can be obtained even in absence of such a detailed knowledge. This 
happens in high energy processes where the only needed nonperturbative input 
are the so-called light-cone wave functions (WF's) with a minimal number of 
constituents. As is known those WF's give parametrically leading contribu- 
tions to hard exclusive processes, hard diffractive electroproduction and many 
others high energy processes. In all such cases the quark and antiquark are pro- 
duced at small distances z ~ 1/Q — ► 0, where Q is the typical large momentum 
transfer. Thus, we can neglect the z 2 dependence of the wave function of the 
meson with momentum p and can concentrate on the variable zp ~ I which is 
order of one. Therefore, the problem is drastically simplified in the asymptotic 
limit and we end up with the light cone wave function which depends only on 
one variable 4>(zp, z 2 = 0). 

The corresponding wave functions have been introduced to the theory in 
the late seventies and early eighties u in ordec. to describe the exclusive pro- 
cesses in QCD. We refer to the review papers0pl,Q on this subject for the detail 
definitions and discussions in the given context. The main idea of the approach 
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is the separation of the large and small distance physics. At small distances 
we can use the standard perturbative expansion due to the asymptotic freedom 
and smallness of the coupling constant. All nontrivial, large distance physics is 
hidden into the nonperturbative WF in this approach. It can not be found by 
perturbative technique, but rather should be extracted from elsewhere. The 
most powerful analytical nonperturbative method for such problems is the 
QCD sum rules a, Q. Therefore, problem becomes tractable within existing 
nonperturbative methods which are based on the so-called Wilson Operator 
Product Expansion (OPE). Therefore, the problem is reduced to the analysis 
of the bound states within OPE. The subject of the Section 2 is the corre- 
sponding analysis of the light cone WF's within QCD. 

We have therefore a well-formulated problem of the high energy behavior. 
The formulation is based on the solid background of QCD. However, calcu- 
lations of the pQCD approach refer, strictly speaking, only to asymptotically 
high energies. Therefore, as usual, the main question remains: at what 
energies do the asymptotic formulae start to work? 

This problem has been a controversial issue for the last fifteen years. In 
order to answer on the question formulated above, we should know the power 
corrections to an amplitude under consideration. It is commonly believed that 
the power corrections are related to the transverse momentum distribution in 
a hadron. Therefore, any dependence on gives some power corrections to 
the leading terms. Naively one may expect that these corrections should be 
small enough already in the few GeV 2 region. However, this is not the case. In 
the Section 3 we develop a theory of the transverse momentum distribution 
in a hadron. This is a key element of the whole approach. The answer on the 
question about k\ will give us the answer on the formulated above question 
regarding the onset of the asymptotic regime. 

The subject of the Section 4 is an application of the these, apparently 
pure academic results to the phenomenological needs of the particle physics. 
Namely, we shall see that the onset of the asymptotic regime is very different 
for different processes in spite of the fact that WF's which come into the game 
are always the same. 

It is instructive to recall a history of this development to the problem 
about an applicability of the QCD-based approach to exclusive amplitudes. 

As we already mentioned, the dependence on k\ gives some power correc- 
tions to the asymptotically leading terms. The expectation that these correc- 
tions should be small enough already in the few GeV 2 region is mainly based 
not on a theoretical analysis but rather on the phenomenological observation 
that the dimensional counting rules, proposed in early seventies, seel3, agree 
well with the experimental data (such as the pion and nucleon form factors, 
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large angle elastic scattering cross sections and so on). This agreement can be 
interpreted as a strong argument that power corrections are small in the few 
GeV 2 region. 

However, in mid eighties the applicability of the approachEJ at experimen- 
tally accessible momentum transfers was questioned u, u. In these papers it 
was demonstrated, that the perturbative, asymptotically leading contribution 
is much smaller than the nonleading "soft" one. Similar conclusion,Jiimporting 
this result, came from the different side, from the QCD sum rules E^j, where 
the direct calculation of the form factor has been presented at /3ir-^ 3GeV 2 . 
This method has been extended later for larger Q 2 < lOGeV 2 Ii3pil with the 
same qualitative result: the soft contribution is more important in this inter- 
mediate region than the leading one. 

Therefore, nowadays it is commonly accepted that the asymptotically lead- 
ing contribution to the exclusive amplitudes cannot provide the experimen- 
tally observable absolute values at accessible momentum transfers. If we go 
along with this proposition, then the natural question arises: How can one 
explain the very good agreement between the experimental data and dimen- 
sional counting rules if the asymptotically leading contribution cannot explain 
the data for experimentally accessible energies? A possible answer was sug- 
gested recently til and can be formulated in the following way: very unusual 
properties of the transverse momentum distribution of a hadron lead to the 
mimicry of the dimensional counting rules by the soft mechanism at the ex- 
tended range of intermediate momentum transfers. Numerically, the soft term 
is still more important than the asymptotically leading contribution at rather 
high Q 2 ~ 50- IQOGeV 2 . 

A situation with hard diffractive electroproduction of the p meson is quite 
different. Having the experience with exclusive processes in mind, one could 
expect a similar behavior (i.e. a very slow approach to the asymptotically 
leading prediction) for the diffractive electroproduction as well. We shall argue 
in the Section 4 that this naive expectation is wrong. The asymptotically 
leading formula starts to work already in the region Q 2 ~ 10 GeV 2 , where 
power corrections do not exceed the 20% level in the amplitude. This surprising 
result is a consequence of very special properties of the WF in the k\ variable 
(as we mentioned earlier, the k\ dependence determines corrections to the 
leading term). 

In Section 5 we discuss and formulate some unsolved problems: How one 
can understand the constituent quark model in terms of QCD (more specifi- 
cally, in terms of .QJJD vacuum structure)? We discuss the so-called Instan- 
ton Liquid Modell 15 !, 16 ! which in its present form correctly reproduces multiple 
mesonic/baryonic/glueball correlation functions, and also has an increasing di- 
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rect support from lattice studies of instantons (see for revie 423). We find that 
many properties of the light cone WF's which have been known for a while, 
can be understood from the instanton point of view. 

2 kj_) and its longitudinal distribution. 

2.1 General constraints 

The aim of this section is to provide necessary definitions and establish some es- 
sentially model independent constraints on the non-perturbative WF kj_) 
which follow from the use of such general methods as dispersion relations, 
duality and large order perturbation theory. 

We define the it meson axial wave function in the following gauge- invariant 
way: 

if7Tq^(zq,z 2 ) = (0|eZ(z)7 M 7 5 e 8fl ^ ^- g Attdz " ' u{-z)\-K{q)) (1) 

= E -MmiM™« A,) n «(0)k(g)>, 
' n! 

n 

where D V =D V — D v and iD^ = id^+gA"^ is the covariant derivative. From 
its definition is clear that the set of different it meson matrix elements defines 
the nonperturbative wave function ij}^{zq, z 2 ). 

The most important part (at asymptotically high q 2 ) is the one related 
to the longitudinal distribution. In this case z 2 ~ the WF depends only on 
one zq- variable. The corresponding Fourier transformed wave function will be 
denoted as <fi(£) and its n— th moment is given by the following local matrix 
element: 

(0| d 7u75 (i D M = iUqv(zq) n (n = %J,q v (zq) n J dK n <t>(0 (2) 

-<7 2 ^oo, zq~l £ = 2u-l, z 2 = 0. 

Therefore, if we knew all matrix elements (^) ( which are well-defined ) we 
could restore the whole distribution amplitude <j>(£). In the infinite momentum 
frame (IMF) q z — * oo the distribution amplitude (DA) </>(£) describes the 
distribution of the total longitudinal momentum q z between the quark and 
antiquark carrying the momenta uq z and (1 — u)q z , respectively. In what 
follows we will use the both variables £ and u , u = 1 — u interchangeously. 

The QCD sum rules approach allows one to find the magnitudes only 
the few first moments 113. As is known, this information is not enough to 
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reconstruct the WF; the parametric behavior at £ — > ±1 is the crucial issue in 
this reconstruction. 

To extract the corresponding information, we use the following duality 
argument. Instead of consideration of the pion DA itself, we study the following 
correlation function with pion quantum numbers: 

i [ dxe i<lx {Q\TJl{x), J o (0)|0) = (zq) n+2 I n (q 2 ), 4 = d lv z vlb {i D M z^) n u (3) 



and calculate its asymptotic behavior at large q 2 . The result can be presented 
in the form of the dispersion integral, whose spectral density is determined by 
the pure perturbative one-loop diagram: 

- [°° ds 11711 "^ , Iml n (sf ert = — %, -. (4) 

7ri s-q 2 ' " W 47r(n + l)(n + 3) 1 ' 

We assume that the tt meson gives a nonzero contribution to the dispersion 
integral for arbitrary n and, in particular, for n — > oo. Formally, it can be 
written in the following way 

1 fS" i roo 

- / dsImI(s)P ert = - / dsImI(s)Z, (5) 
tt Jo ' t Jo 

Our assumption means that there are no special cancellations and it meson 
contribution to the dispersion integral is not zero, i.e. S™(\\) ^ 0, where we 
specified the notation for the longitudinal distribution. In this case at q 2 — > oo 
our assumption (|^) leads to the following relation: 

£<rx»->oo)->^ (6) 

It unambiguously implies the following behavior at the end-point region i: 

(H = £ dt?<KZ) ~ l/n 2 , ±1)-(1 -?)■ (7) 

Few comments are in order. Because we consider (by definition) the total 
correlation function (0), the behavior (j|) should be fulfilled for any nonper- 
turbative WF no matter what the specific shape of WF is. Thus, our first 
constraint looks as follows: 

.1 0(£->±l)^(l-£ 2 ). (8) 
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We want to emphasize that the constraint ( »1) is of very general origin and 
follows directly from QCD. No numerical approximations were involved in 
the above derivation. Pre-asymptotic as q 2 — > — oo perturbative and non- 
perturbative corrections are only able to change the duality interval in Eq. (||) 
(which is an irrelevant issue, anyhow) but not the parametric 1/n 2 behavior 
which remains unaffected. 



2.2 Numerical constraints in the longitudinal direction 

We shall discuss here the 7r meson case only. Similar consideration can be 
carried out far other hadrons as well within the same technique. We refer to 
review articlaa for the references references. 

The QCD sum rules for the moments of the tt meson WF defined above 
(0) have the following form: 

i A" dse-^ 2 lml n (s) = ff ' + (9) 

7r J A-K z {n + l)(n + 3) 

1 ,a ari2 . 167r(ll + 4n) 2 



'-GU+ 81M4 H V^qqY 



-Iml n (s) = f 2 (C)S(s) + 6(s-s n ) 



4ir 2 (n + l)(ra + 3)' 



where M 2 is so called Borel parameter which varies in the region where power 
corrections as well as continuum contribution (modeled by the standard 9(s — 
s n ) function) are small enough, i.e. < 30%. The result of the standard fitting- 
procedure with respect to the Borel parameter M 2 looks as follows: 



(eV = 1.5GeV- = J i dti 2 <f>(0 - 0.38, (£V=2.2GeV- ~ 0.21 



(10) 



Now, in order to interpret these results, we should compare them with the 
corresponding values of the asymptotic DA (f>(^) asym = 3/4(1 — £ 2 ): 

($ 2 ) = J d£€ 2 4>(t)' uym - 0.20, (e 4 ) = 0..086 (11) 

For comparison of those DA's we should keep in mind that: 

a) They both normalized to unity, d£(j)(£,) asyrn = j\ d£</>(£) = 1, 

b) They both have the same behavior at ^ — > ±1, (f>{^ — > ±1) — > (1 — ^ 2 ), 

c) They have very different few first moments, see eqs (|Io|), (|lT|). 

Such a result was a motivation to suggest the so-called " two-hump" or "CZ" 
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wave functional. Very soon, such a \iLF -became a very popular and controver- 
sial issue at the same time, see, e.g. Ej"o. 

It is not the purpose of these lectures to comment the corresponding re- 
sults. Rather, I would like to make a remark, that there are strong arguments 
in favor as well as against of this two-hump WF. My opinion is: the true re- 
sult rests somewhere in the middle of those two numbers: © and dn]). In 
different words, the QCD sum rule predictions (|Io|) gi ve somewhat larger mag- 
nitudes for the moments and asymptotic results (11) predict somewhat lower 
values in comparison with what we believe would be the correct magnitude of 
the moments. 

Unfortunately the information regarding (£ 2 ) is very difficult to extract by 
doing a phenomenological analysis of an experiment. The problem in such of 
analysis is the following. Before to extract any information from the experi- 
ment, we have to make sure that one of the following conditions is satisfied: 

a) . We are already in the asymptotic regime, i.e. the power corrections are 
small. Therefore, an asymptotic formula already works and one can extract a 
relevant information from the corresponding theoretical expression. 

b) . We are not in the asymptotic regime. However, the power corrections are 
under control, i.e. we know the origin of those power corrections as well as we 
know a way of how to estimate them. In this case, again, one could extract 
the relevant information regarding the longitudinal distribution </>(£). 

Unfortunately, none of those conditions is satisfied presently in any exper- 
iment. Therefore, one can not make a reliable prediction about the properties 
of the distribution amplitude 0(£) unless we know precisely the power correc- 
tions. That is why an interpretation of an experiment could be very ambiguous 
and very often it is based on a strong assumption about power corrections. 

We can not improve our understanding by going along the condition a). 
To satisfy this condition we need to go to the much higher energy and momen- 
tum transfer in comparison to what is available at the moment. Therefore, we 
shall go along the line b), i.e. we shall try to understand the power correc- 
tions relevant to the problem. As we mentioned in Introduction, those power 
corrections are tightly connected to the behavior of a k\) in transverse 
direction k\. Actually it was the main motivation for thfi, recent study of the 
hadron transverse momentum distribution within QCEcJp^3. We present a 
review of those results in the next section. 



3 Transverse distribution 
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3.1 General constraints 



The moments in transverse directions are defined analogously to Eq. (g) through 
gauge invariant matrix elements 



C2r? - IV - 

(o\d 7 , 15 (t d v t„y n u\K( q )) = iM-t 2 r [ (2n)! /" (fcf) , (12) 

where transverse vector = (0, t, 0) is perpendicular to the hadron momentum 
Qfj, = {QOi 0j_, q z ). The factor ^g^fr is introduced to ( |l2| ) to take into account 
the integration over </> angle in the transverse plane: J d(f)(cos 4>) 2n / J d<fi = 
(2n — l)!!/(2n)!L By analogy with a non-gauge theory we call (kj_) in this 
equation the mean value of the quark perpendicular momentum, though it 
does not have a two-particle interpretation. Indeed, it is very different from 
the naive, gauge dependent definition like (0\dj l/ j5d 2 ^u\Tr(q)} , because the phys- 
ical transverse gluon is a participant of this definition. We believe that such 
definition is the useful generalization of the transverse momentum conception 
for the interactive quark system. Its relation to the higher Fock components 
will be discussed at the end of this section. Its relation to constituent quark 
model and Instanton picture of QCD vacuum will be discussed in Section 5. 
Here we note that Eq. ( |l2| ) is the only possible way to define (kj_) in a manner 
consistent with gauge invariance and operator product expansion. 

To find the behavior (fcj™) at large n we can repeat the previous duality 
arguments with the following result []: 

fW± )^^ - " ! =* f*$±) - - ! (13) 

This behavior has been obtained in RefS by studying of the large order per- 
turbative series for a proper correlation function. Dispersion relations and 



duality arguments transform this information into Eq.(13). It is important to 
stress that any nonperturbative wave function should respect Eq.(|l3|) in spite 
of the fact that apparently we calculate only the perturbative part ( see the 
comment after Eq. ([?])). The duality turn this perturbative information into 
exact properties of the non-perturbative WF. The most essential feature of 
Eq. ( |l3| ) is its finiteness for arbitrary n. This means that higher moments 

(kl n )= I dk\d^{kl,0 (14) 



Here and in what follows we ignore any mild (non-factorial) n-dependence. 
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do exist for any n. In this formula we introduced the non-perturbative 
kj_) normalized to one. Its moments are determined by the local ma- 
trix elements ( |l2| ) which are obviously finite, because they are normalized at 
low normalization point /i. The relations to Brodsky and Lepage notations 
^ bl{x\, fcjjtj and to the longitudinal distribution amplitude (f>(£) introduced 
earlier, look as follow: 

9 BL (xt,k x ) = ^^*(S,ki), J dkl^,kl) = m, J^iO = 1 

(15) 

where f n — l33MeV. The existence of the arbitrary high moments (fc^™) 
means that the non-perturbative ^5>(£,k^), defined above, falls off at large 



transverse momentum k\ faster than any power function. The relation (13) 
fixes the asymptotic behavior of ^(£,, k±) at large k\. Thus, we arrive to the 
following constraint: 



• 2 (kT) = J dkid^kf 1 ^^, k 2 ± ) ~ n! n -> oo. (16) 

We can now repeat our duality arguments again for an arbitrary number of 
transverse derivatives and large (n — > oo) number of longitudinal derivatives. 
The result reads S 



.3 y"dfcifci fc *(fci,e^±i)-(i-e 2 ) fc+1 (it) 



The constraint (|17]) is extremely important and implies that the k]_ dependence 
of *if(kj_,£) comes exclusively in the combination kj_/(l — £ 2 ) at £ — > 
±1. This means that the standard assumption on factorizability \&(fc^,£) = 
xp(kj_)(fi(£) is at variance with very general properties of the theory such as 
duality and dispersion relations. The only form of kj_) satisfying all the 
constraints (|7]),(|l6|) and ( |l7| ) is the Gaussian with a very particular argument: 

vK(fci - oo,f - ±1) ~ exp ^ A2(1 fc ^ 2) ^) (18) 

(here A 2 is a mass scale which can be fixed by calculating the moments 
(kj_),(k±) etc.) Strictly speaking, so far we have only established the va- 
lidity of Eq.(nq) in a vicinity of the end-point region £ — > ±1. However, one 
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can argue that , the behavior (18) can be approximately valid (with some 
accuracy) in the whole range of the £ variable. 

We would like to pause here in order to make the following conjecture. 
The Gaussian ^(£, kj_) (reconstructed above from the QCD analysis) not acci- 
dentally coincides with the harmonic oscillator iff(£,kj_) from the constituent 
quark model. We shall discuss this conjecture within the instanton model in a 
more detail in Section5, but now let us recall some results from the constituent 
quark model. 

It has been known for a whilc3 that the equal- time (ET) wave functions 
*bt(9 2 ) ~ exp(-<f 2 ) (19) 

of the harmonic oscillator in the rest frame give a very reasonable description 
of static meson properties. Together with Brodsky-Huang-Lepage prescrip- 
tion E2l connecting the equal -time and the light-cone wave functions of two 
constituents (with mass m ~ SQOMeV) by identification 



one can reproduce the Gaussian behavior (18) found from QCD. It means, first 
of all, that our identification of the moments ( Q ) defined in QCD with the 
ones defined in quark model, is the reasonable conjecture. The same method 
can be applied for the analysis of the asymptotic behavior of the nucleon WF 
which in obvious notations takes the form: 

£2 

IWZeon^L -> OO) %i) ~ 6Xp(- ^ — )■ (20) 

However, there is a difference. In quark model we do have a parameter which 
describes the mass of constituent m ~ 300Afey. We have nothing like that 
in QCD. This difference has very important phenomenological consequences 
which will be discussed in the next section. Here we would like to emphasize 
that there is no room for such mass term in QCD. Its inclusion violates the 
duality constraint ^ since in this case we would have 

(D ~ I <%e exp (- ™* „\ ~ exp(-v^) , n ^ co (21) 
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A 2 (i-e 2 



instead of the 1/n 2 behavior (^). In other words, a true nonperturbative WF 
must respect the asymptotic freedom which is incompatible with the quark 
model -type mass term in the WF. To set this more accurately, one can say 
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that a possible (scale-dependent) mass term in *f>(£, k±) must renormalize to 
zero at a normalization point p 2 ~ a few GeV 2 where the duality arguments 
apply. This conclusion is not at variance with popular models for the QCD 
vacuum such as e.g. the instanton vacuum (see til for review). Even more 
than that: The instanton picture does support all properties we have been 
discussing in these lectures, see the section 5 for details. In no sense we claim 
that a WF like exhausts the tt ( p) -meson properties, or (j|^) exhausts 
the nucleon properties. 

3.2 Lowest moments (kj_) , (kj_) an d vacuum condensates. 

The general constraints of the previous section are insufficient for building up 
a realistic non-perturbative WF for the p, n -mesons or nucleon. To fix this, 
we follow the same logic as in the analysis of the distribution amplitudes and 
calculate few lowest (fej™},n = l,n = 2 moments of k\). The physical 
meaning of the second moment (kj_) is clear: this quantity serves as a common 
scale for power corrections in physical amplitudes. The fourth moment (kj_) 
tell us on how strongly kj_) fluctuates in the k\ plane. In this section we 
will calculate the lowest moments (kj_) and (k\) for p and tt meson kj_). 
With the knowledge of and (kj_) we then construct a model WF , J , (£, k 2 ) 
which will be used to estimate higher twist effects in different processes in the 
next section. 

We start our discussion from the the second moment of the kj_) in the 
transverse direction-defined by equation ( |l2] ) . For the tt meson it was calculated 
for the first time in El and independently (using a quite different technique and 
very different motivation) inH3. Both results are in a full agreement to each 
other: 

/f2 \ ^ ^9 a p.vG°^ v \q) 5to§ 2 2 2 / 00 n 

{k±)n = — 7^7 ^ -^r ^ (330Mel/) , m ~ 0.8GeV . (22) 

A very similar calculation can be done for the p meson as welEI. The result 
looks like this: 

(k 2 ± ) p = (420MeV0 2 (23) 



Essentially, the result given by equations (|22|) and (|23|) defines the general scale 
of all nonperturbative phenomena for the pion and p meson correspondingly. 
Those numbers are very much the same. It is not accidentally coincides with 
300 -7- 400MeV scale which is the typical magnitude in the hadronic physics. 

To study the fine properties of the transverse distribution it is desired to 
know the next moment. The problem can be reduced to the analysis of the 
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mixed vacuum condensates of dimension seven 

8 4(gq) 9(gq) 

A similar formula for the p meson has been derived in0 and takes the following 
form: 

/p4\ 3 4/ 4\ , 1 (0|gg 2 G Mt /G^g|0) 1 (Ol^cr^G^axaGx^qlQ} 
{k ± ) p _ -m p <« ) + j g (25) 

Therefore, we have explicitly expressed (kj_) in terms of the vacuum expec- 
tation values (VEV's) of the dimension 7 operators. This is very important 
result: we essentially say that the properties of a hadron are determined by 
the QCD vacuum structure! One can say even stronger: the gluon fluctuations 
inside of a hadron are essentially nothing, but the vacuum fluctuations. They 
are very strong (they have in general lGeV scale) but they are approximately 
the same for all hadrons. This observation gives a chance to understand a 
constituent quark model: all hadrons are build up from the same constituents 
which originated from QCD vacuum, see section 5 for details. 

Coming back to (]24|) and (|H) one could estimate these condensates naively, 
by factorizing them into the products of the quark (^V) an d gluon (g 2 G 2 ) con- 
densates. This procedure, based on the factorization hypothesis, does not work 
in the given case: there arc essential deviations from the factorization predic- 
tion in VEV's (^3[|4|). The non-factoraziblity, of mixed quark-gluon matrix 
elements of such type has been studied incilH by two independent methods 
with full agreement in estimates between them. The first onecil was based 
on the analysis of heavy-light quark systems, while the second method has 
related the vacuum condensates of the form (|2^^) to some pion matrix ele- 
ments known from PCAC. Here we only formulate the result of this analysis. A 
measure of non-factorizability is introduced by the correction factors K\ , K 2 
in the matrix elements 



1 

A 

G 

■1 



WG^G^O) = -^i(0|«7 2 G^G^|0) (0\qq\> 



(0\qg 2 a^G^a Xa G Xa q\0) = —K 2 (0\g 2 G 2 \0) (0\qq\0) (26) 

( K\ — K2 = 1 in the factorization limit). Those factors are approximately the 
same K\ = K2 = K ~ 3 for the both mixed operators appealing in Eq. (p6|) . 
A possible uncertainty of this estimate does not exceed 30 %E3. In this case, 
formulae (p5j), ( p6[ ) give the following estimate: 

(ki) p ~ O.UGeV 4 . (27) 
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We would like to present the result of these calculation for (kj_) p by introducing 
a dimensionless parameter which is quantitative measure of the fluctuations in 
the transverse direction: 

R= 0^~ 4 ~5. (28) 

For the tt and p mesons parameter R is almost the same and it is very large. 
It corresponds to an unexpectedly large hadronic matrix elements of the high 
dimensional operators. We have explicitly calculated those m.e. for the oper- 
ators of dimension seven. However, we expect that a similar conclusion also 
takes place for arbitrary higher dimensional operators.. Such a large value of 
the parameter R corresponds to the strong fluctuations in the transverse di- 
rection. In terms of wave function this property means a very unhomogeneous 
distribution in the transverse direction. 

In terms of the QCD vacuum structure these fluctuations are due to the nu- 
merical enhancement of the high dimensional quark gluon mixed condensates 
or, what is the same, the large magnitude of the parameter nonfactorizability 
K in Eq. ( |2"rj| ) . On the microscopical level such an unhomogeneous distribution 
corresponds to some fluctuations of the strong gluon fields with a small size. In 
this case, factorization prescription (which essentially suggests a homogeneous 
distribution of the gluon vacuum fields) clearly does not work. As a conse- 
quence of this, high dimensional operators do violate a factorization. 
We identify these vacuum fluctuations with instantonsEZI, see section 5. 

An additional evidence in favor of this view-point is the recent calculation 
of the rf matrix element (which can be measured experimentally!) of a high- 
dimensional gluon operator within an instanton model. In general, it is very 
difficult to have a direct experimental measurement of a matrix element of a 
high dimensional operator. In most cases, we use some indirect experiments 
(see section 4: Applications) to establish the properties of those matrix ele- 
ments. Fortunately, due to the uniqueness of the rf meson a high dimensional 
operator can be measured in the CLEO directly! Such a matrix element can 
be explicitly extracted from the B — > rf decays^. This matrix element 
indeed is very large in agreement with our analysis. 

As we shall discuss in section 5, in the constituent quark model those 
properties correspond to the small size (lGeV -1 ) of the constituent quark in 
comparison with the hadron size (1/m). 

3.3 Model wave function k\) 

The results obtained so far are essentially model independent. We have fixed 
the form of the high-£r[ tail of the true nonperturbative WF k\) Eq.([l8|) 
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by the use of the quark-hadron duality and dispersion relations. Furthermore, 
we calculated the lowest moments (fcj_) and (k]_) using the equations of motion 
and QCD sum rules. Now our purpose is to build some model for the true 
non-perturbative WF fcjj which would respect all general constraints of 
the previous Sections and incorporate the effect of strong fluctuations in the 
transverse k\ plane found in the previous section (see Eq.(|2g|)). 

We start our discussiprLfroni the analysis of the \I/(£, kj_) motivated by 
constituent quark modelEj,E]~LJ (CQM) Q Such a function is known to give a 
reasonable description of static hadron properties. The Brodsky-Huang-Lepage 
prescriptions leads to the following form for the pion ^(^, k 2 ^): 

r*2 2 

*(kl,u) C QM=Aexp(-^±^), (29) 

We call this function as the constituent quark model WF. As we already dis- 
cussed before, it satisfies two general constraints (•2,»3), but not to (»1) be- 
cause of the nonzero magnitude for the constituent mass m. We make the stan- 
dard choice for the parameter m ~ 330MeV in accordance with its physical 
meaning. The parameter /3 can be determined from the numerical constraint 



for the mean value (kj) for ir meson (22) and p meson (|23| ) correspondingly. 



They do not differ much and we shall not distinguish them for qualitative 



discussions. Parameter A is determined by the normalization eq.(15) 



To make function wider in the longitudinal direction, one can insert into 



formula (29) an additional factor 



(l + <7(M)K 2 -k .9(m)=5(M2)-(^4) 50/9b (30) 

with additional parameter g(p). This new parameter g(p) allows to adjust (£ 2 ) 
as appropriate. For tt and p meson WF's this parameter is quite different. For 
the asymptotic distribution amplitude parameter g = 0. 

Now we go from Constituent Quark Model to QCD-based WF. Before to 
design QCD{k\, u), let us explain what do we mean by that. We define the 
non-perturbative wave function ^(fc^, u, p)qcd through its moments which can 
be expressed in terms of the non-perturbative matrix elements ([l)Jl2]). As is 
known, all non-perturbative matrix elements are defined in such a way that 
all gluon's and quark's virtualities smaller than some parameter \i (point of 
normalization) are hidden in the definition of the " non-perturbative matrix 

c Here we neglect all terms in QCM related to spin part of constituents. In particular, we 
do not consider Melosh transformation and other ingredients of the light cone^equal time 
connection. It does not effect any qualitative results presented here. 
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elements" (|l2]). All virtualities larger than that should be taken into account 
explicitly (perturbatively, due to the asymptotic freedom). In particular, all 
perturbative tails like l/kj_ should be subtracted from the non-perturbative 
WF by definition. The same procedure should be applied for the calculation 
of non-perturbative vacuum condensate {G J), where the perturbative part 
related to free gluon propagator 1/fc 2 is also should be subtracted. 

With these general remarks in mind we propose the following form for the 
non-perturbative wave function ^{k 2 ^, u, ^o)qcd at the lowest normalization 
point: 

y(ki,u^o) Q CD = Ae M-^^) ■ {1 + .9(Mo)[e 2 - \}}. (31) 

Parameter (3 for this parametrization is found to be j3 ~ 0.3 GeV (it corre- 
sponds to R = 2.2 and (fc^) = O.lAGeV 2 ). In comparison with the constituent 
quark model the "only" difference is the absence of the mass term ~ m in the 
exponent. This difference is a key element. We discussed this point earlier and 
we would like to emphasize this point again: the nonzero mass in ^(fc^, u)cqm 
was unavoidable part of the wave function within a quark model. As we dis- 
cussed earlier we do not see any room for such a term in QCD, because its 
presence corresponds to the behavior 

<o = J ^em ~ J der* ex p(-Tz^ ) ~ e M-M, n - «>, (32) 

which is in contradiction to 1/n 2 result found earlier (Q). 

One can implement the effect of strong k\ fluctuations discussed previously 



and which is quantitatively expressed in terms of parameter i?(28). It can be 
done in a number_of ways. We shall not discuss this point here refering to the 
original literatures. However, we note, that this effect leads to the existence of 
two characteristic scales. in kj_) which has its explanation within instanton 
vacuum liquid model 113. As we already mentioned, such a property of the 
could be viewed as an explicit manifestation of the complexity of 
the hadronic structure: one scale ~ Ifm determines the hadron size itself; 
another, new scale ~ \GeV~ x determines the hadron substructure, the size of 
the constituent quark. 

Our conclusion is that the transverse momentum distributions for the 7r 
and p-mesons are to a large extent alike. However, the distribution in the 
longitudinal direction is very different for those WF's. 

Let us summarize. We constructed two different types of wave functions. 
The first one, ^&cqm is motivated by quark model with its specific mass pa- 
rameters. The second type is motivated by QCD consideration. All these 
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WFs have Gaussian behavior at large k\. However, in the case of ^cqm this 
behavior is related to the nonrelativistic oscillator model, while for QCD moti- 
vated models this behavior is provided by constraints discussed in the previous 
section. 

Contrary to the CQM, the QCD motivated wave functions do not contain 
the mass parameter m ~ 300Af eV which is an essential ingredient of any quark 
model. Such a term is absolutely forbidden from the QCD point of view. 

In the next section we discuss some applications. In particular, we calcu- 
late the contribution to the pion form factor caused by these wave functions. 
We shall find the qualitative difference in behavior on Q 2 , which is our main 
point. We shall also discuss the behavior of the hard diffractive electroproduc- 
tion of p meson as a function of Q 2 . Again, as in it meson form factor case we 
shall see a strong influence of the transverse momentum distribution on this 
behavior. We shall also discuss nucleon form factor. 

4 Applications 

4-1 7r meson form factor 

The first application we considerjs the pion form factor. The starting point 
is the famous Drell-Yan formula Cj (for modern, QCD- motivated employing 
of this formula, seeB), where the F^(Q 2 ) is expressed in terms of full wave 
functions: 

/dxd 2 k i 
-i^r^Bifo k± + (1 - x)$l)*bl(x, k ± ), (33) 

where q 2 = — q\ = —Q 2 is the momentum transfer. In this formula, the 
1 $>BL{x,k±) is the full wave function; the perturbative tail of &bl(x, k±) be- 
haves as a s /fc^ for large k\ and should be taken into account explicitly in 
the calculations. This gives the one-gluon-exchange (asymptotically leading) 
contribution to the 7r meson form factor in terms of distribution amplitude 

In terms of QCD, the formula ( |33| ) is an assumption. However, we ex- 
pect, that by taking into account only "soft" gluon contribution (hidden in 
the definition of k\ (El)), we catch the main effect of the soft physics. There 
is no proof for that within QCD. The only argumentation which can be deliv- 
ered now to support this assumption is based on the intuitive picture of quark 
model, where the current quark and soft gluons form a constituent quark with 
original quantum numbers. No evidence for the gluon playing the role of a 
valence participant with a finite amount of momentum is seen. 
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In terms of vacuum structure those soft gluons are nothing but vacuum 
fluctuations. Therefore, those gluons are very strong in amplitude (they are 
soft in a sense of the momentum they carry on which is small). Those vac- 
uum fluctuations (which are classical configurations like instantons) are much 
stronger than the quantum fluctuations carrying a nonzero momentum. 

From the viewpoint of the operator product expansion, the assumption for- 
mulated above, corresponds to the summing up a subset of higher-dimension 
power corrections. This subset actually is formed from the infinite number of 
soft gluons and unambiguously singled out by the definition of non-perturbative 

We refer to the original paperli-l for details. Here we formulate the result 
of calculations. The main qualitative difference between quark model and 
QCD- motivated wave functions is as follows: a much slower fall off at large 
Q 2 is observed for the k\) motivated by QCD. The qualitative reason for 
that is the absence of the mass term, see discussion after the formula (|32|). 
Precisely this term was responsible for the very steep behavior in all previous 
calculations based on a quark model wave function. The declining of the 
form factor getting even slower if one takes into account the property of the 
broadening of ^{^,k\) in transverse direction. This property corresponds to 
the strong fluctuations in the transverse direction and quantitatively is related 
to the large parameter R discussed in the previous section. 

Therefore, our main observation is that the QCD based WF's could mimicry 
the dimensional counting rules by the soft mechanism at the extended range 
of intermediate momentum transfers. Numerically, the soft term is still more 
important than the asymptotically leading contribution at rather high Q 2 ~ 
50 WOGeV 2 . 

4-2 Nucleon form factor 

Now we would like to extend our previous analysis to the nucleon form factor. 
The starting point, as before, is the fundamental constraints (»1 — »3), which 
being applied to the nucleon wave function imply the Gaussian behavior with 
the specific argument (^p|). With these constraints in mind one can model 
the nucleon WF in the same way as we did for the pion. Having modeled 
the nucleon wave function, one can calculate a soft contribution to different 
nucleon amplitudes. The corresponding analysis was carried out in the ref. E3. 
Here we quote some results from this paper. 

The most important qualitative result of these calculations is similar to 
what we already observed previously in the 7r- meson case: namely, the com- 
bination Q 4 F nucl (Q 2 ) is almost constant in the extent region of Q 2 in spite 
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of the fact that the corresponding "soft" contribution naively should be de- 
creasing function of Q 2 . The qualitative explanation of this phenomenon is 
the same as before and is related to the absence of the mass term in the QCD- 
motivated wave function. 

The next observation is related to the longitudinal distribution and can be 
formulated as follows: A fit to the different experimental data leads to a wave 
function which has the same type of asymmetry which was found previously 
from the QCD sum rules. The asymmetry is however much more moderate 
numerically than QCD sum rules indicate. 

In particular, one can calculate the valence quark distribution functions 
u p (x) and d p (x) at large x in terms of the non-perturbative nucleon WfE3. 
Two properties of the WF are important to provide such a fit: The absence 
of the mass term in the formula ( p0| ) (this leads to the correct power behavior 
at x — > 1) and a moderate asymmetry in the longitudinal direction ( this 
provides an observed ratio for the ^ p |^j ~ 5 at x — ► 1 in the contrast with the 
asymptotic formula prediction which gives value of 2 for the same ratio). 

We already mentioned earlier that the similar conclusion is likely to have 
place for the pion wave function also. Therefore, the general moral, based on 
already completed calculations, can be formulated in the following way. There 
is a standard viewpoint for the phenomenological success of the dimensional 
counting rules: it is based on the prejudice that the leading twist contribution 
plays the main role in most cases. This outlook, as we mentioned earlier, is 
based on the experimental data, where the dimensional counting rules work 
very well. We suggest here some different explanation for this phenomenolog- 
ical success. Our explanation of the slow falling off of the soft contribution 
with energy is due to the specific properties of non-perturbative kj_). In 
particular, we argued that the absence of the of the mass parameter in the 
corresponding kj_) is the strict QCD constraint. At the same time this 
property is responsible for the behavior mentioned above. Besides that, we 
found a new scale (~ lGeV 2 ) in the problem, in addition to the standard 
low energy parameter (k\) ~ O.lGeV 2 . Both these phenomena lead to the 
temporarily mimicry of the leading twist behavior in the extent region of Q 2 . 

We believe that this is a new explanation of the phenomenological success 
of the dimensional counting rules at available, very modest energies. Besides 
that, as we shall see in the last section, the new scale we found has very jiatiixal 
explanation in terms of the specific Vacuum configurations, instantonffll 
They have small size ~ 1/3 fm which is very different from the hadron size 
~ 1/m. The last scale is determined by the density of instantons. The new 
1/3 fm scale, we believe, should appear in the quark model. However, in 
the constituent quark model this new scale can be nothing, but the size of a 
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constituent quark, see section 5 for details. 



4-3 Hard diffractive electroproduction 

In this section we apply our model WF for the study of the pre-asymptotic 
effect due to the Fermi motion in diffractive electroproduction of the p-meson. 
Our prime goal is to get an estimate for the onset of the asymptotic regime in 
this problem. 

The applicability of perturbative QCD (pQCD) to the asymptotic limit p£ 
the hard diffractive electroproduction of vector meson was established in RefLj 
using the light-cone perturbation theory. The authors of Ej have proved that 
for the production of longitudinally polarized vector mesons by longitudinally 
polarized virtual photons the cross section can be consistently calculated in 
pQCD. There was found that at high Q 2 the amplitude factorizes in a product 
of DA's of the vector meson and virtual photon, the light-cone gluon distri- 
bution function of a target, and a perturbatively calculable on-shell scattering 
amplitude of a qq pair off the gluon field of the target. After the factorization 
of gluons from the target the problem seems to be tractable within OPE-like 
methods. 

Let us start with the general form of the amplitude as a matrix element 
of the electromagnetic current = e(2/3u7 M u — l/Sdj^d) : 

M = e^N(p - r)p(q + r)\j^q)\N(p)) , (34) 

where e M is the polarization vector of the photon (only the longitudinal polar- 
ization is considered, seeEj) and r stands for the momentum transfer. We will 
neglect the masses of the nucleon and p-meson in comparison with the photon 
virtuality Q 2 : m 2 p , m 2 N <C Q 2 . For the momentum transfer r we consider 
the limit r 2 — 0, but r M ^ 0. Wc jiext note the following. The factorization 
of the gluon field of the target EJ, Ell means that that these gluons act as an 
external field on highly virtual quarks produced by the photon with Q 2 — ► oo. 
Retaining only the leading contribution, we arrive to the asymptotic formula 



of RefEl in the form suggested in 
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4nV2ea s f p 1 f 1 lv F°(X)VT=I #«) 



N c QJ X(X-Z + i5)J uu V ' 

where (j>(u) is the standard p meson light-cone DA and F^(X) is so-called asym- 
metric gluon distribution which becomes the usual gluon distribution function 
XFg(X) in the symmetric limit £ — > 0. 

The most important corrections to the asymptotic formula (^Hj) are due 
to the quark transverse degrees of freedom (the Fermi motion). Calculating 
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only this contribution we make an educated guess on the scale of higher twist 
corrections in the diffractive electtoproduction. Technically, this correction can 
be written in the following formE^j 

fl du cQ 2 mZ y( p \ 1 ( i _ o fcj/(ml) \ 

, JO UU JO " KJ -*V "> «-jJ (Q 2+fc2 /(ufi)) 2 1 1 Z (Q2 + fc2 /(u0) I 

v / W) = <9 4 ; ^ — ^ .> 

(36) 

where ^(u, fc^) was defined earlier in terms of the local matrix elements ( j bi| ). 
By definition, for the asymptotically large Q 2 — > oo we have \/T(Q 2 — » oo) = 
1. Deviations from \/T(Q 2 ) = 1 determine a region of applicability of the 
asymptotic formula (pq). 

Now we are in position to discuss numerical estimates for the correction 
factor (|3^) in order to answer the main question formulated above. The authors 
offi3 have observed that the choice of ^(w, fejj in a factorized form ^(u, kj_) — 
(j>(u)ip(k±) leads to a very slowly raising function \/T(Q 2 ) which approaches 
0.8 at rather high Q 2 > 40 GeV 2 depending on the model chosen for ^(u, kj_). 
Using the QCD motivated WF we see a very fast approach to the asymptotic 



in (£56|). The correction factor \/T(Q 2 ) reaches the value 0.8 already at Q 2 ~ 
10 GeV 2 . Such a behavior implies a rather low onset for the asymptotic 
regime for the hard diffractive electroproduction. The technical reason for this 
behavior is quite clear: large values of k\/(uu) are exponentially suppressed 



by the QCD motivated WF and thus these extra terms in thti integral (36) 
give a small contribution for Q 2 > 10 GeV 2 , see original pape£J for details. 

Thus our final conclusion is that the onset of the asymptotic regime for 
diffractive electroproduction of the longitudinally polarized p-meson is approx- 
imately Q 2 ~ 10 GeV 2 where corrections due to the quark transverse degree of 
freedom constitute less than 20 % in the amplitude. This can be traced back to 
the fact that in the case at hand the power corrections are given by the matrix 
elements of local operators and in fact are fixed completely by the independent 
calculation of the moments. This is the consequence of the structure of Eq.(^6|) 
and the fact that the WF is a function of the single variable kj_/(uu). Only 
global, but not local characteristics of ^(u, k\) in the k]_ plane are important. 
This situation can be confronted with the case of exclusive processes. In the 
most well studied problem of the pion form factor, the , asymptotic regime has 
been found to be pushed further to Q 2 3> 10 GeF 2 where the sub-leading 
"soft" contribution is still larger than the leading asymptotic one. There are 



d We follow the notations of Refl 53 ! by reserving the symbol T(Q 2 ) for the correction in the 
cross section. 
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no reasons to expect the onsets of the asymptotic regime to be similar in the 
diffractive electroproduction and exclusive processes. On contrary, they differ 
parametrically in l/a s . Moreover, the explicit calculations suggest that the 
asymptotic regime in the p-meson diffractive electroproduction starts already 
at Q 2 ~ 10 GeV 2 . We stress that this conclusion refers only to the diffractive 
electroproduction of the longitudinally polarized p-meson, the situation with 
the transverse polarization or diffractive charmonium production can be quite 
different. 



4-4 Concluding remarks 

To conclude this section we would like to make a comment regarding the shape 
of the 7r meson distribution amplitude. This issue has been discussing for quite 
a while. We already mentioned at the end of section 2 that a direct extraction 
of the corresponding information from the experimental data is an extremely 
difficult (if possible at all) problem. Only indirect analysis is available at the 
moment. The point is that the power corrections (which are very difficult to 
estimate) could be very important at present energies. To be more specific, 
and in order to clarify our point, lei ua consider amplitude 77* — > n which 



has received some attention recentl^' l2EjL3. Some theoretical calculations 



claim that the "two-hamp" wave function disagrees with CLEO results; other 
calculation^ conclude that either the asymptotic or the "two-hamp" wap: 
function is sufficient to describe the data. Much more sophisticated methocEJ 
which avoids any assumptions regardine_a shape of the wave function, gives a 
very good description of the experimenlEfl. ILowever, in spite of the power and 
generality of the method advocated in refEa, it does not allow us to extract 
an information regarding the shape of the wave function. In order to obtain 
such an information we have to make some additional strong assumptions 
regarding the power corrections (as well as perturbative coirections) within 
the standard scheme of PQCD calculations originated in refstl. Therefore, our 
point is that a direct analysis of the available experimental data can not provide 
an unambiguous information about the properties of the wave functions. 



5 Instantons and the Constituent Quark Model. 

Let us remind that the purpose of these lectures is twofold. First, these lectures 
have strong applied (or phenomenological) direction. The main result obtained 
in this direction could be formulated in the following way: there are many theo- 
retical, phenomenological and experimental evidences that the "soft" contribu- 
tion (being un-leading parametrically) nevertheless can temporarily mimic the 
leading twist behavior in the extent region of Q 2 : 3GeV 2 < Q 2 < AOGeV 2 . 
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This is due to the specific properties of ^{kj_,x) we have been discussing at 
length in these lectures. Such a mechanism, if it is correct, would be an expla- 
nation of the phenomenological success of the dimensional counting rules at 
available, very modest energies for many different processes. At the same time, 
the same properties of ^>(kj_, x) applied to the hard electroproduction predict 
the onset of the asymptotic regime at the relatively small Q 2 ~ lOGeV 2 . We 
believe that these are very new results and very new understanding which 
deserve for the further studying. 

The second goal is much more theoretical at the present level of under- 
standing, but it may become a profound phenomenological tool in the hadronic 
physics in future (at least we hope so). We have in mind the,|QCI)-based de- 
scription of the WF properties within the instanton modeEali3tZl. A reader 
probably has noticed that almost in each section^ while we discuss a result 
obtained in a phenomenological way, we also try to explain the same result 
in terms of the QCD vacuum structure which could provide such a behavior. 
We also have been trying to explain those results in terms of the specific (but 
important) vacuum configurations which could be responsible for the aforemen- 
tioned properties. We believe that those vacuum configurations are instantons. 

Now we would like to collect all relevant information from the previous 
sections here, in one place, in order to formulate our understanding of the quark 
model in terms of the instanton vacuum configurations. This presentation 
unavoidably can not be self-contained, but rather it is a very fragmentary one. 

5.1 Instantons. Historical Remarks. 

The instanton solution as the classical solution with nontrivial topology was in- 
vented in 1975 motivated by the discovery of the 't Hooft-Polyakov monopolc. 
Shortly after its discovery, the physical meaning of the instanton as a tunneling 
event between-degenerate classical vacua was understood. We refer to a peda- 
gogical revie\Cj for details. Here we shortly introduce some relevant notations 
by emphasizing on the physical meaning of the instanton solution. 

The starting point is the transition from Minkowski to Euclidean space. In 
this case the best tunneling path in gauge theory which connects topologically 
different classical vacua is the solution of classical equation of motion in Eu- 
clidean space. To find these solutions, it is convenient to exploit the following 
identity 

^ = 4^2 J d^xG^G^v = ^2 J d^x^G^G^ + -(G^ v =p G<* u ) 2 ], (37) 

e Sfie, in particular, a discussion after eq. (|2l[); the end of section 3.2; a discussion after eq. 
(|32[); the end of section 4.2. 
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where G^ v = l/2e ljLl/ \ l jG ( ^ a . Since the first term is a topological charge Q of 
the configuration and must be an integer, se&B 

one can argue that the action is minimal if the field is self-dual or antiself-dual 
= iG" . We call these solutions as an instanton or antiinstanton corre- 
spondingly. From (|37]) we have S — 87r g j,^ for (anti)self-dual fields, implying 

that the tunneling probability is Ptunneiing ~ exp (— 87r ffi ), where the coef- 
ficient in front of the exponent is determined by a one-loop calculation. The 
corresponding calculations have been completed by 't Hooft in 1976. We do 
not need in our discussions an explicit formula of the 't Hooft calculations, 
however we do need to know few general results which follow from his formula, 
sc(E3 for details: 

1. An each instanton is accompanied by 4iV c zero modes (N c = 3 is a 
number of colors), which are related to the so-called collective coordinates 
describing the field of the instanton. Those degrees of freedom are: center of 
the instanton Xq (4 modes) size of the instanton p{l mode) orientations of the 
instanton in the color space f2$, (4JV C — 5 modes). 

2. Density of instantons dnj, or what is the same, tunneling probability is 
divergent at large p: 

P~^ ~ex P (-^)%^(A P )», 6 = lliV c /3 = ll. (39) 
9 KP) P b P b 

Such a result is a consequence of the unjustified perturbative calculations of 
the function at large size p. We come back to this point later. 

3. An each instanton is accompanied by the fermionic zero mode as a 
consequence of the very general so-called Index theorem. 

5.2 Chiral symmetry breaking and instantons. 

It is believed that the one of the most profound features of QCD, the chiral 
symmetry breaking phenomenon, can be explained by instantons. This idea 
has been under discussion for a_jquite a while. The modern approach to the 
problem has been suggested inEj, see also reviewaJ. We refer to the original 
papers on the subject for details, here we formulate the main ideas and results 
of this study. 

As we mentioned above, in the presence of an instanton, there is an exactly 
one zero mode. This zero mode is the chiral mode: it is right-handed for 
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instantons and left-handed for antiinstantons. As we shall see, this fact has 
the profound phcnomenological consequences for many aspects of the theory: 
a description of the chiral symmetry breaking phenomenon, the formation of 
the chiral condensate, there appearance of a nonzero constituent quark mass 
and many others. 

Suppose we have infinitely separated / (instanton) and / (antiinstanton) 
which can be thought as two degenerate states with two zero modes. Like in 
the standard quantum mechanical calculations, when we decrease a distance 
between / and /, this degeneracy is lifted through the diagonalization of the 
hamiltonian. As the result, the obtained two new states have non-zero eigen- 
values X± which are equal to the overlap integral Tjj between the original 
states A± = ±|T 7 j|. When one adds more I's and I's, each of them brings in 
a would be zero mode. After the diagonalization they get split symmetrically 
with respect to zero point A = 0. Eventually, for an instanton- antiinstanton 
ensemble one gets a continuous band spectrum with a spectral density v{X) 
which is finite at A = 0. Such a spectrum implies a nonzero magnitude for the 
chiral condensate (qq) ~ 7w(A = 0)E3. 

Having explained the physical mechanism of the chiral symmetry breaking 
phenomenon we proceed with the discussion of some quantities which are much 
closer to the phcnomenological needs of the hadronic physics. Specifically, 
let us consider the quark propagator in the instanton vacuoim we discussed 
above. The corresponding formula has been derived in refO. The result of 
this calculation has the form of a massive propagator with a momentum- 
dependent dynamical mass: 

s{k) = p + wS) ' M{k) = MoJr2(fc) ' M °- 350MeV ( 4 °) 

F{z)=2z[I (z)K 1 {z)-I 1 (z)K {z)--I 1 (z)K 1 {z)], F(0) = 1 (41) 

z 

In this formula z = l/2kp (p is the size of instanton, p — 1/3/to); M(k — > 
oo) ~ -p^- — * as it should due to the asymptotic freedom. However, at small 
k, the mass is not zero but finite number M(k — * 0) = Mo ~ 350MeV which is 
close to the phenomenological value for the constituent quark mass. The value 
of the chiral condensate is also can be calculated in this model in a similar way 
with the following result 

which is again very close to the phenomenological magnitude. We conclude this 
short historical introduction into instanton physics with the following general 
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remark: A systematic numerical study of various correlation functions, ampli- 
tudes, hadronic characteristics (like charge radius of a hadron, constant / w , etc) 
in the instanton vacuum demonstrates a good agreement with experiments and 
phenomenology, see review£3. Similar conclusion has been recently obtained 
from direct lattice measurements. Therefore, it is fair to say that instan- 
tons explain the basic properties of light hadrons and large distance dynamics. 
With this optimistic note about instantons we proceed with the discussion of 
our main subject: the light cone wave functions. 



5.3 Low energy: how to formulate the problem? 

We start from the following simple observation: At low energies it is very dif- 
ficult to define a WF in terms of the original quark and gluon fields. The 
problem is that infinite number of different matrix elements are equally impor- 
tant in such a would be definition. All of them are not zero and all of them 
have a standard hadronic scale IGeV. Therefore, there is no any special rea- 
son to select a quark- aniquark WF in comparison, let us say, a quark-aniquark 
and ten -gluon -field WF. To be more specific and in order to explain the main 
point, let us consider few possible choices. One could start from the simplest 
two-particle WF (we consider pseudoscalar WF for simplification) which is 
normalized to the following matrix element: 

, . — . . fn'm'i (uu) + (old) 

di l5 u 7T = J * * c~ - V ' V 1 . (43) 
m u + m d }n 

One could make one step further and add an arbitrary number of gluons into 
this definition: 



A\„ . , 2(g(i g <7 M „G%„^)"g) 

) U ITT > = 

2 ' ' ' h 



(0|&7S(i5^G>— ) n «k) = ~ — t ^ 2J H/ , (44) 



where in the last equation we have used PCAC to reduce the matrix element 
to the vacuum condensates. All those vacuum condensates are not known 
exactly, but it is believed that they are not zero and have a normal hadronic 
scale which is about IGeV. One could define a WF, based on these matrix 
elements ( ff4| ) in the same way as we did in eq.(Q): non of them is better or 
worse than others at low energy. Moral of this exercise is simple: There is no 
any reason to start from the simplest WF at low energy, because an arbitrary 
number of gluons are not suppressed at all. Therefore, all of them go on the 
same footing. 
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5.4 High energy: twist classification 



Situation at high energies, as we explained in Introduction is quite different: 
there is a twist classification which makes possible to select the leading twist 
WF uniquely. Specifically, for the ir meson WF it is given by eq.(P with 
the axial Lorentz structure. All other WFs give parametrically smaller con- 
tribution ~ (l/Q 2 )™ to any amplitude. Therefore, we can neglect all of them 
except the leading ones, which by dimensional reasons, are exactly the WFs 
with a minimal number of constituents. In different words, precisely those 
WFs represent a constituent nature of the hadrons. 

What happens when we go to the lower and lower energies? We do not have 
the complete answer on this question, but we hope we convinced a reader (see 
section 5.1) that the most important vacuum gluon fluctuations in hadronic 
low-energy physics are the instantons. Once we accept this, we know the 
answer on the question formulated above: The most important quantitative 
change which instantons make is the complete reconstruction of the relevant 
degrees of freedom: gluons have disappeared; instead a constituent quark with 
a momentum-dependent mass appears M(k) = AfoF 2 (fc)(ffO|). 

5.5 Next steps. 

Having explained the non-perturbative physics which is related to the instan- 
tons, we now in position to answer on the question (at least qualitatively) 
formulated above: what happens to the leading twist WFs when energy is 
getting smaller and smaller? Our proposal for the answer is: All those WFs 
acquire a momentum-dependent mass term. In particular, ~^{k\, u, IIq)qcd 
( |3l| ) takes the form: 

^{kl lU ,^) QCD = Aexp(- x W ) • {1 + gMie - ■=]}, (45) 

with function M(z) given by eq.([fO|). A similar replacement should be made 
for all WFs we have discussed previously. The main feature of this replacement 
is clear: At high energies when a typical k is the same order of magnitude as 
the external Q, we have M{k\ — + oo) — ► 0. In this region, WF ( |45|) transforms 
into the leading twist WF given by formulae (|3l|). At small energies when 
M(kj_) — > Mo, our WF ( (45|) transforms into the constituent quark model WF 
given by eq.(p8|). Most important question: What happens to a variety of WFs 
we mentioned above (H) in such a transition? Where do they go? 

Before to answer on this question, we should remark that a large magnitude 
for all hadronic matrix elements with gluons (similar to (0)), is related to the 



2G 



strong vacuum fields (which we identify specifically with instantons in this 
section) and not to a specific properties of a hadron. In the case of the 7r 
meson this relation follows from the PCAC. For different hadrons a similar 
relation is not so obvious as for pion. Nevertheless, experience with p meson 
showSZl, that all gluonic matrix elements are large and they are related to the 
vacuum condensates, i.e. to the vacuum fields. If it is so, all would be the new 
WFs (similar to eq. (ff4|)) describe nothing, but some vacuuum fluctuations. 
The main effect related to those vacuum fluctuations is well known: a quark 
becomes a constituent quark with the momentum dependent mass M(k). In 
different words, vacuum gluons are transformed into different powers of k. 

Therefore, the answer on the formulated above question regarding the 
numerous of WFs is: they transform into the only relevant degrees of freedom, 
the massive strongly interacting constituent quarks with dependent on k mass. 
In terms of these constituent quarks, an apparent variety of those WFs simply 
disappears. This effect can be explained using an analogy with the harmonic 
oscillator problem in the external electric field. As is known, in this case, the 
whole problem is reduced to the change of variable (one should make a shift 
of the coordinate on the amount which is proportional to the external electric 
field). It is exactly what happens in our case: variable k\ get shifted on the 
amount of M(k 2 1 _). Instanton fields are not constant fields, therefore, the value 
of a shift ~ M depends on momentum k\ . This analogy is even closer to the 
real situation, because a Gaussian dependence which we derived from the QCD 
earlier (see section 3.1) is exactly the solution of the Schrodinger equation for 
the harmonic oscillator potential. 

Now we can reverse our arguments in order to explain the reason for the 



matrix elements (44) to be large. Those gluon fields G^ v which are present 
in the definition of the hadronic matrix elements have nothing to do with the 
specific hadronic state, but rather, they are related to some strong vacuum 
(instanton) fluctuations. In this picture we also understand some relations, 
like (Pq ) which have been obtained in a pure phenomenological way. Namely, 
the fluctuations in the transverse direction for the different hadrons alike f\ and 
numerically those fluctuations are very large. In the picture we suggest this 
phenomenon has its natural explanation: All matrix elements which describe 
transverse momentum distribution are related to the vacuum (instanton) fields; 
therefore they are the same for any hadron. Secondly, large magnitude for the 



parameter R(28) is related to the nonfactorizability of the vacuum condensates 
which is also has a natural explanation from the instanton point of view. New 
scale we discussed in the previous sections is nothing but the size of the in- 



•^Wc checked this explicitly only for the tt and p mesons, but we believe that a situation will 
be the same for all hadrons. 
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stanton p ~ 1/3 /m in the formula for M (fc)([Io|). This scale is much smaller 
than the quark mass Mq ~ 350MeV. 

To end this section let me remark that WF (^) with the momentum- 
dependent mass may help to resolve many long-standing problems in the in- 
termediate region of the nuclear and low-energy particle physics. With this 
hope we conclude our lectures. 
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